Two-scale analysis of the SU(N) Kondo Model 
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We show how to resolve coherent low-energy features embedded in a broad high-energy background 
by use of a fully self-consistent calculation for composite particle operators. The method generalizes 
the formulation of Roth, which linearizes the dynamics of composite operators at any energy scale. 
Self-consistent equations are derived and analyzed in the case of the single-impurity SU (N) Kondo 
model. 



71.10.-w, 71.27.+a, 75.20.Hr 

The development of effective methods for describing 
correlated electron systems has been the subject of in- 
tensive activity over the last decade, spurred by the ex- 
perimental discoveries of the heavy fermion systems, the 
high-temperature superconductivity and, generally, a re- 
vival of interest in transition metal-oxide physics jjj . 

The Roth method for the correlation problem Q| in the 
context of the Hubbard model j|] is based on an ansatz 
which reduces the dynamics of field operators to a lin- 
earized one. The essential idea is to select a basis of 
fermionic operators ipi, write their equations of motion 
which involve operators Ji, and then close these equa- 
tions by projecting Ji onto the basis by using the Roth 
projector V defined by 

V{J l )=Y,({Jl^l})lrsHs (1) 

rs 

where I rs = ({ipr, V 7 ]})! with {_,_} denoting the anti- 
commutator. In this approach the determination of the 
Green's functions is then reduced to the evaluation of cer- 
tain static thermal averages: ({yj r , V'l}) an d {{Jh V4})- 
When these parameters are connected to matrix elements 
of the Green's functions associated to the basis, one has 
a self-consistent scheme for their calculation. However, 
this is often not the case, and further approximations are 
introduced for their evaluation. The application of this 
method, as well as similar methods [Q, has recently been 
reviewed by Mancini and collaborators Q . Through care- 
ful comparison with existing numerical data, they con- 
cluded that good results for many physical quantities are 
obtained by requiring that the Green's functions fullfil 
exact equal-time identities accompanying the fermionic 
character of the operators. 

In spite of its intuitive appeal, there are several serious 
difficulties with the Roth's method. Recent advances in 
the study of correlated electron systems converge upon 
a picture of the one-particle Green's function made up 
of incoherent broad spectral features in addition to more 
dispersive quasi-particle bands which exist at lower ener- 
gies ||,|7j . The Roth approach describes the Green's func- 
tions in terms of a finite number of sharp poles which are 
a poor description of the incoherent structure of the high- 
energy spectra. Also, the presence of low-energy features 



embedded in a broad high-energy background precludes 
the straightforward extension to low-energy scales of this 
approach. Indeed, low-energy features cannot be resolved 
increasing the size of the basis. Increasing the size of 
the basis only amounts to calculating self-consistently a 
larger number of spectral moments |8| which are domi- 
nated by high-energy contributions. A clear example of 
this dramatic failure is provided by the Kondo impurity 
model, where it has been proved impossible to derive the 
existence of a Kondo resonance in the spectra within a 
projection scheme. 

In this Letter, we present a generalization of Roth's 
projection technique which overcomes the limitations dis- 
cussed above and, as an illustration of the technique, we 
investigate the single-impurity SU(N) Kondo model 
Our goal is to introduce a general technique in a simple 
context which is well understood, but so far has not been 
successfully treated by techniques based on the equations 
of motion. We will demonstrate that our method repro- 
duces all the well-known spectral features of the impurity 
model. 

The method carries out the following steps, (i) In the 
first step, we write the equations of motion for the opera- 
tors of physical interest in terms of higher order ones (or 
composite operators). Similarly, we express the Green's 
functions of interest in terms of the Green's functions 
of the composite operators. The composite operators 
should not have components on the physical fields at 
high energies, (ii) Then, we evaluate the Green's func- 
tions of the composite operators by a technique which 
is valid at high energies, such as the mode-coupling ap- 
proximation Use of the mode-coupling approxima- 
tion is motivated by the fact that it more clearly reflects 
that the high-energy part of the spectra is quite incoher- 
ent. At this stage, we divide the composite operators into 
(a) a high-energy part, described by the mode-coupling 
approximation, and (b) a low-energy part, to be deter- 
mined by a non-perturbative closure of the equations of 
motion. The necessity for this step is checked by writ- 
ing the expressions for the moments and noting that the 
mode-coupling approximation fails to give the spectral 
weights, as expected from an independent evaluation of 
the moments, (iii) The low-energy closure of the equa- 
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tions of motion is dictated by the physics of the problem 
and is inspired by the successes of the slave boson tech- 
niques |ll]] . It is a simple quasi-particle theory which in- 
volves unknown parameters such as the low-energy spec- 
tral weights. The self-consistent determination of the 
low-energy parameters completes the full determination 
of the physical Green's functions. 

The SU(N) Kondo model is described by the following 
Hamiltonian: 



ff = XV(k) 
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where c(k) denotes the conduction electron operator, and 
n d represents the spin operator at the impurity site (n d ■ 
n d = N(N + l)/2). e c (k) and Jk are the conduction 
electron energy and the Kondo coupling, respectively. N s 
is the number of atomic sites of the host metal responsible 
for the orbitals which form the conduction band, rfc are 
the iV 2 — 1 traceless generators of SU(N) (tn ■ fjv = 
2(N 2 -l)/N). From (|), we have 



i| C (k)=e c (k)c(k) + 2J K ^ 
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Next, we introduce the Composite Heisenberg Field 
Operator: 
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where Co = S q c ( t l) i s the electron at the impurity 
site. The field 02 m (H) abides by the criterion required 
by the method in (i). In fact, when ^|'02, - 0ij - ^ is re- 
garded as the scalar product of the field 02 with the field 
Co, there is no component of 02 on cq at any energy scale 
since ^{02,0ij^ = 0- Then, using the equation (||) 
we can express the Green's functions of the first field in 
terms of the Green's function for the composite operator 
■02- We have 

G u ( w ) = roH + 4r HG 22 Hr H (5) 

G 12 (oj)= J K r (w)G 22 (a/) (6) 

where G a p{oj) is the thermal Green's function associated 
with the basis in (^) and Tq(uj) — ^_ 1 \ \ is the 
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free propagator of the field Co. The total spectral weight 
attached to the second composite field 02 is 



I22 = 4 



N + 1 

N 2 



(7) 



where the Kondo amplitude Kb = ^0i0|^ = 

— 2jj? (ip\ ■ tn 0i n d ^j describes the binding between the 
localized spin and the spin excitations of the field cq . In 



order to resolve low-energy features embedded in a high- 
energy background, we write 



G22M = Gf 2 H + G&M 



(8) 



where G^w) keeps the information about the band 
structure and is not sensitive to features which are small 
with respect to the bandwidth 2D. In contrast, G^w) 
mostly takes coherent contribution from low energies and 
depends only weakly on the high-energy part of the spec- 
trum. Such a decomposition corresponds to decomposi- 
tion of the composite field 2 as 02 = 1P2 + ^2, with 
1P2 giving rise to incoherent broad features, whereas ip% 
emerges as an observable quasi-particle at low energies. 

In the high-energy regime, time-dependent correlation 
functions can be treated within the mode-coupling ap- 
proximation in terms of electron-hole and charge-spin 
fluctuations. By use of mode-coupling in the paramag- 
netic case, we have for the time ordered Green's function 
52 H: 
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where 

Sii r (ti,tj) = {T [n d r {U) nf(tj)] ) r = 1, . . . , N 2 - 1 (10) 

The spectral weight absorbed by the propagator G22<w) 
in the mode-coupling form (|J) is 4(iV + 1)/N 2 . For sim- 
plicity, we take the atomic limit for the Bose propagator, 
so that G 2 c 2 (w) = 4-^trGii(cj). Other treatments for 
the Bose propagators would not substantially affect our 
results for the fcrmionic spectral function. 

From the Hamiltonian (|^) it is direct to derive 
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where f^ bc are the structure constants of the SU(N) 



Lie algebra ([t#,t^] 
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Tf^Eq £ c(q)c(q). It is worth noting that the source 
(O) has a direct component on the field cq (2Jk/Ntn • 
^rip2 = 4Jk(N + 1)/N 2 cq + .. .) which disappears for 
N -> 00 being the coefficient AJ K (N+ l)/N 2 (for N = 2, 
2J K /Nf N -n d ^2^ 3J K c - 2J K 2 ). 

In the low-energy regime, we assume the following dy- 
namics for the field ip2- 



0_ 
'df 
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This corresponds to the physical assumption that at low 
energies (i.e., at energies much smaller than Jk and D) 
we have a quasi-particle theory. Indeed, the ansatz in 
( |l2| ) can be described as an application of Roth's pro- 
jection idea to a field 02 which has most of its spectral 
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weight at low energies. Thus, the high-energy spectral 
weight is already accounted for by the mode-coupling ap- 
proximation. This is the second main departure from the 
original Roth approach, where an equation of motion for 
the field ip2 would have been projected onto ip\ and 1^2 it- 
self. This field would have spectral weight at all frequen- 
cies (e.g., 4(N + 1)/N 2 ) and from it a Kondo scale cannot 
be estimated. Once again, noteworthy is the fact that the 
basis defined in ^ alone is inadequate to capture both 
low- and high- energy physics of the Kondo model once 
a Roth truncation is realized. At this level of approxi- 
mation, in the scattering matrix (i.e., J^G 22 (lo) ) there 
is only one energy scale that cannot mimic a crossover 
between the two regimes. This is set by I22 where the 
high-energy spectral weight (i.e., 4(iV + 1)/N 2 ) prevents 
the low-energy scale from emerging. 

By combining (j|) and ([l2|) it is direct to show that 
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being K = J^I 22 , after projecting (|l2| ) on the field ip\. 
We have defined I 22 as the spectral weight of G 22 {oj) in 
the low-energy region. Also, it is implicitly assumed that 

^|V , 2,'02 /t }) = because they span different energy 
sectors of the Hilbert space. In conclusion, we have 
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At this stage of the method, once I 22 i s evaluated, the 
problem is solved, as in point (hi) referred to above. From 
Eq. (Q) it is clear that I 22 is connected to Kjj- However, 
it is crucial to note that it represents the low-energy spec- 
tral weight and is a distinct object in respect to the Kjj 
parameter which also contains contributions from ener- 
gies of the order Jk- For the determination of this low- 
energy scale, we need to call upon some self-consistent 
condition. Evaluating the Kondo amplitude Kjj using 
G12M 



K D = -Ty^G 12 (f u> n ) 



(2n + l)nT (17) 




FIG. 1. The Kondo amplitude Kn is plotted as a function 
of the temperature for Jk = 0.08 and N — 2. 




FIG. 2. The low-energy spectral weight I22 is given as a 
function of the temperature for Jk = 0.08 and N = 2 (solid 
line). As a point of reference, we give the solution after 
Jg = F [I22 = 0] (dashed line) and for iV = 00 (dotted line). 



and inserting this into Eq. (Q), we give a relation be- 
tween ^22 and I 22 of the form ^22 = F [l 22 ] ■ To 
get the self-consistent equation, we estimate I 22 = 
F [l 22 = 0, T = Jk] which results in an equation for the 
low-energy spectral weight of the form I 22 = ^22 — I22 • in 
other words, we choose Jk as the energy above which the 
one-particle Green's function is made up of incoherent 
high-energy contributions with no relevant temperature 
dependence. While the equations resemble the slave bo- 
son equations, they differ in a significant way from them. 
These equations do not introduce additional redundant 
phases and Lagrange multipliers, and avoid all the diffi- 
culties associated with the treatment of the gauge fields. 
The slave boson method has been very successful in ob- 
taining low-energy information. High-energy information 
can also be obtained by performing fluctuations around 
the mean-field solutions, but this gets increasingly diffi- 
cult particularly in lattice models |12] . 

We now present some results. For the numerical 
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FIG. 3. The spectral density a-niyi) is reported for 
T = 0.001 (solid line) and T = 0.1 (dashed line). J K = 0.08 
and N = 2. 

solution of the equations we used a constant density 
of states p = jp8 (D — \u\) for the field Co so that 
r«( W ) = ^ln|g+H|-i^CD-M)- In Figs. [J and 
0, the Kondo amplitude Krj and are shown as func- 
tions of the temperature for Jk = 0.08 and N = 2. 
D has been set equal to 1. In Fig. || we also show 
the solution for N = oo and the one after replacing 
I22 = F [J& = 0,T = J K ] with J* = F [J& = 0] . Both 
these solutions give a spurious transition at a character- 
istic temperature which is the signature of the Kondo 
crossover in the exact solution. As in the slave boson 
approximation, this is due to the absence of a small in- 
homogeneous term which is present in our method and 
mimics mixing effects between high- and low- energy con- 
tributions. The quantities in Fig. || coincide only in the 
limit of large spin degeneracy (i.e., N — > oo) where the 
method recovers the exact results H . Our numerical esti- 
mate for the Kondo temperature Tk agrees well with the 
exact solution |p|,^3[. At this temperature (Tk — 0.002), 
I22 has a change in the concavity of its slope when plotted 
as a function of the temperature. In Fig. 0, we present 
the spectral density er 2 2(w) = (— — ) 3 [G^w)] for two 
different temperatures. Again, it is clear that at high 
temperatures (i.e., T > Jk) only a high-energy incoher- 
ent background is left. A well-defined singlet excitation 
mode no longer exists, even if some residual spin-spin in- 
teraction persists being energetically favored by a finite 
bandwidth (i.e., Kd is non-zero at any temperature as 
in Fig. 0). 

In conclusion, we have shown how to resolve coher- 
ent low-energy features embedded in a broad high-energy 
background by use of a fully self-consistent calculation for 
composite particle operators. In a problem with more 
than one energy scale, which is typical of strongly corre- 
lated systems, we succeeded to capture low-energy fea- 
tures. Our scheme extends and improves upon the Roth's 
method by combining the advantages of the methods 



based on the equations of motion and the slave boson 
techniques. Finally, we note that when there is an expan- 
sion parameter such as the size of the group, or the size 
of the representation, our approach can be formulated so 
as to reduce to the correct solution in the exactly solu- 
ble limit. We have illustrated this here with the Kondo 
model in the limit of large spin symmetry group which 
has often been shown to retain many crucial aspects of 
low-energy physics Q . 

Finally, our approach is directly applicable to lattice 
models and work in this direction is currently in progress. 
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